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SUMMARY 

A unified  prediction  method  has  been  developed  to  support  V/STOL 
Stability  and  Control  analyses.  The  method  is  geared  to  a preliminary 
design  environment  to  facilitate  its  application  using  the  limited  data 
which  are  normally  available  early  in  the  aircraft  design  process. 

Mathematical  models  have  been  formulated  for  the  force  and  moment 
contributions  of  each  of  the  vehicle  components  and  then  combined  to  allow 
for  prediction  of:  non-linear,  six  DOF  (Degrees  of  Freedom)  trim  conditions; 

three  DOF,  small  perturbation  stability  characteristics;  six  DOF  response 
time  histories  for  arbitrary  control  and  disturbance  inputs.  The  structure 
of  the  methodology  follows,  in  a gross  sense,  a similar  methodology  for 
helicopter  and  stoppable  rotor  aircraft  (reference  (a)),  with  considerable 
modification  of  the  individual  component  models.  Primary  emphasis  has  been 
placed  on  the  jet  lift  and  lift/cruise  fan  classes  of  V/STOL  vehicles  in 
the  hover  and  transition  flight  regimes  but  the  methods  should  be  easily 
extendable  to  include  other  classes  of  vehicles. 

The  methods  have  been  programmed  for  the  CDC  6600  computer  and  have 
been  validated,  with  reasonable  success,  using  available  XV-6A  Kestrel 
and  AV-8A  Harrier  static  and  dynamic  data.  Results  of  the  validation 
effort  are  presented  in  Appendix  A. 

Documentation  is  in  two  volumes.  This  volume,  Volume  I,  describes 
the  component  models,  analytical  methods  and  validation  results. 

Volume  II  constitutes  a User  Manual  for  the  digital  computer  program. 
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INTRODUCTION 

The  stability  and  control  characteristics  of  V/STOL  aircraft  in  hover 
and  transition  flight  can  be  quite  different  than  those  of  CTOL  aircraft  in 
the  corresponding  take-off  and  power  approach  flight  regimes.  V/STOL  aircraft 
are  usually  unstable  at  low  speeds  and  experience  propulsion  induced 
aerodynamic  and  inlet  momentum  drag  effects  which  are  normally  neglected 
for  CTOL  aircraft.  The  problem  is  further  complicated  by  additional  control 
functions  such  as  primary  engine  nozzle  angle  and  differential  thrust  and 
RCS  (Reaction  Control  System)  thrust  modulation.  A recent  survey  of  V/STOL 
aircraft  designers  (reference  (b))  has  indicated  a serious  need  for  V/STOL 
stability  and  control  prediction  methods  to  support  preliminary  design, 
configuration  trade-offs  and  evaluations.  To  this  end,  the  development 
described  herein  was  undertaken. 

Many  past  efforts  in  this  field  have  been  limited  to  uncoupled, 
two  or  three  DOF  analyses  at  small  angles  of  attack  and  sideslip.  Since 
V/STOL  vehicles  may  exhibit  significant  cross-coupling  between  longitudinal 
and  lateral  axes,  it  is  necessary  to  consider  full  six  DOF  simulation  with 
cross-axis  coupling.  Angles  of  attack  and  sideslip  from  0°  to  +180°  may 
be  encountered  in  low-speed  flight  requiring  complete  non-linear  aerodynamic 
modelling.  Additionally,  provision  must  be  made  for  adequate  modelling  of 
aerodynamic  interference  effects  resulting  from  propulsion  system  efflux 
interaction  with  the  surrounding  vehicle  surfaces.  All  of  the  above 
effects  and  other  less  significant  ones  were  included  in  the  model  formulation. 
Figure  1 summarizes  the  major  model  components. 

The  prediction  technique  has  been  developed  primarily  for  application 
in  a preliminary  design  environment  and,  hence,  does  not  require  extensive 
wind  tunnel  or  propulsion  system  data  tables.  This,  of  course,  has  some 
impact  on  the  accuracy  of  the  results  but  validation  with  flight  test  data 
for  an  existing  vehicle  (XV-6A  Kestrel)  has  exhibited  an  adequate  degree 
of  accuracy  for  preliminary  vehicle  design,  configuration  trade-off  and 
design  evaluation  from  a stability  and  control  viewpoint.  The  methodology, 
in  its  current  status,  is  applicable  to  lift  plus  lift/cruise  and  lift/cruise 
turbojet  and  fan  configurations.  However,  the  general  model  structure  is 
such  that  extension  to  include  additional  vehicle  classes  is  certainly 
feasible. 

The  aerodynamic  model  of  the  airframe  is  divided  into  four  major 
sections:  fuselage,  wing,  horizontal  tail  (conventional  aft  tail  or  canard) 

and  vertical  tail.  Each  component  is  modelled  separately  using  a combination 
of  theoretical  computations  and  empirical  curve  fits  over  the  complete  angle 
of  attack  and  sideslip  range.  Aerodynamic  effects  which  are  not  directly 
related  to  one  of  these  sections  (e.g.,  landing  gear  effects)  are  included  in 
the  fuselage  simulation  for  simplicity.  Propulsion  induced  aerodynamic 
interference  effects  are  accounted  for  by  using  curve  fits  of  estimated  and 
available  experimental  data. 

A maximum  of  two  fixed  jet  nozzles  and  six  vectorable  jet  nozzles 
may  be  simulated,  with  control  being  maintained  over  the  thrust  of  all 
Jets  and  nozzle  angles.  In  addition  to  the  main  jet  nozzles,  a maximum  of 
ten  RCS  nozzles  and  their  associated  control  functions  may  also  be  simulated. 
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FIGURE  1:  VEHICLE  MODEL  COMPONENTS 
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All  component  models  are  combined  to  form  the  total  vehicle  model  which 
may  then  be  analyzed.  Analytical  procedures  have  been  developed  for  the 
following : 

1.  Non-linear,  six  DOF  trim  solution; 

2.  Six  DOF  stability  derivative  and  control  power  estimation; 

3.  Small  perturbation  characteristic  root  and  transfer  function 
calculations; 

4.  Non-linear,  six  DOF  time  history  response  to  control  inputs 
and  disturbances; 

5.  Least  squares  time  history  parameter  vector  analysis  (amplitude 
and  phase) . 

Details  of  each  component  model  and  analytical  procedures  are  described 
in  the  following  sections.  Discussion  of  input  data  development  is 
presented,  where  necessary,  to  aid  in  the  proper  application  of  the 
prediction  methods. 

The  total  method  has  been  programmed  in  FORTRAN  for  the  CDC  6600 
computer.  Input  and  output  descriptions  and  guidance  for  program  application 
are  included  in  Volume  II. 


MODEL  DESCRIPTION 

The  total  aircraft  model  consists  of  an  integrated  combination  of 
four  major  submodels.  Each  of  the  submodels  and  their  associated  components 
are  listed  below. 

1.  Aerodynamic  Characteristics 

a.  Fuselage  aerodynamics 

b.  Lifting  surface  aerodynamics  (wing,  horizontal  stabilizer, 
vertical  stabilizer) 

c.  High-lift  flap  aerodynamics 

d.  Wing  downwash 

e.  Propulsion-induced  aerodynamics 

2.  Aircraft  Mass  and  Inertia  Characteristics 

3.  Propulsion  System  Characteristics 

a.  Fixed  nozzle  engine  characteristics 

b.  Vectorable  nozzle  engine  characteristics 

c.  Inlet  momentum  effects 

d.  Angular  momentum  effects 
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4.  Control  System  Characteristics 

a.  Aerodynamic  control  surfaces 

b.  Reaction  control  nozzle  characteristics 

c.  Thrust  management  control  characteristics  (magnitude  and 
direction) 

Forces  and  moments  produced  by  each  component  are  calculated  in  a 
localized  axis  system  and  then  resolved  to  a fixed  set  of  reference  aircraft 
body  axes  with  the  X-axis  directed  forward,  Y-axis  directed  toward  the 
right  wing  tip  and  Z-axis  directed  down.  The  centers  of  action  of  each 
component  force  and  moment  system  are  located  in  the  aircraft  reference 
system  by  FS  (Fuselage  Station)  increasing  in  the  negative  X direction, 

BL  (Butt  Line)  increasing  in  the  positive  Y direction  and  WL  (Water  Line) 
increasing  in  the  negative  Z direction.  Each  component  force  and  moment 
system  is  transferred  to  the  aircraft  center  of  gravity  and  all  are  summed 
to  yield  the  total  force  and  moment  acting  on  the  vehicle. 

Each  component  model  will  now  be  described  in  detail. 

AERODYNAMIC  CHARACTERISTICS 
Fuselage  Aerodynamics 

Aerodynamic  characteristics  of  the  fuselage  are  modelled  in  a body 
axis  frame  of  reference  as  axial,  side  and  normal  force;  pitch  and  yaw 
moments.  Roll  moment  is  not  an  explicit  model  input  but  evolves  from 
applying  side  force  through  the  vertical  moment  arm  between  the  fuselage 
aerodynamic  reference  center  and  the  vehicle  center  of  gravity.  Figure  2 
shows  the  orientation  and  sign  convention  of  the  fuselage  forces  and 
moments . 

Since  one  of  the  objectives  of  the  prediction  method  is  to  simulate 
vehicle  characteristics  throughout  hover  and  transition,  force  and  moment 
models  for  the  complete  angle  of  attack  (a)  and  sideslip  (ay)  range  are 

required.  In  the  case  of  the  fuselage,  it  is  assumed  to  be  sufficient  to 
concentrate  on  the  0 to  +180  degree  range  with  the  following  assumptions 
applied  to  account  for  the  0 to  -180  degree  range. 


A (-a)  - A (a)  (la) 

S (-ciy)  = -S  (cty)  (lb) 

Nf  (-a)  --Nf  (a+2a0)  (lc) 

M (-a)  --M  (a+2a0)  (Id) 

N (-aty)  - -N  (Oy)  (le) 


Inherent  in  equations  (1)  are  the  assumptions  that  side  force,  S,  and  yaw 
moment,  N,  are  zero  at  ay  *■  0,  axial  force,  A,  has  a minimum  at  a = = 0, 
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FIGURE  2:  FUSELAGE  FORCE  AND  MOMENT  CONVENTIONS 
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and  normal  force,  Nf,  and  pitch  moment,  M,  are  zero  at  a = a0.  The 
aerodynamic  angles  a and  ay  are  defined  by: 


a 


a 


y 


-1 

tan 


-1 

tan 


iw/u) 

(2a) 

(v/u) 

(2b) 

where  u,  v,  and  w are  the  body  axis  velocities  in  the  X,  Y,  and  Z directions, 
respectively. 

The  aerodynamic  forces  and  moments  of  the  fuselage  must  be  determined 
from  experimental  data  or  estimated  by  theoretical  or  semi-empirical  means. 
Reference  (c)  describes  a very  convenient  digital  computer  technique  for 
estimating  the  aerodynamic  characteristics  of  symmetrical  bodies  over  the 
entire  angle  range.  It  also  includes  a very  good  bibl iography  of  other 
currently  available  techniques.  Figure  3 presents  the  results  of  applying 
the  method  of  reference  (c)  to  the  XV-6A  (Kestrel)  fuselage.  These  data 
were  used  during  validation  which  is  described  in  a later  section.  The 
data  are  normalized  by  dynamic  pressure,  qQ,  to  allow  for  varying  flight 
velocities. 


The  data  curves  of  Figure  3 are  typical  of  fuselage  aerodynamics  and 
may  be  reasonably  approximated  by  the  equations  presented  in  Table  I. 

Though  the  equations  may  seem  somewhat  complex,  a closer  inspection  reveals 
that  determining  the  required  coefficients  is  really  quite  simple. 

(A/q0)0  is  the  value  of  the  axial  force  at  zero  angle  of  attack.  ( )max 

is  the  maximum  value  of  the  corresponding  force  between  0 and  180  degrees. 

( )maxj  an^  ( )raax2  are  t*ie  maxinuim  values  of  the  corresponding  moment 


between  0 and  90  degrees  and  90  and  180  degrees,  respectively.  Note  that 

(M/q0)  and  (N/q0)  are  usually  negative.  The  moments  take  on  zero 

° max2  max2 

value  at  some  angle  between  0 and  180  degrees;  these  angles  are  denoted  by 

for  pitch  and  <*2  for  yaw.  Once  the  preceeding  coefficients  are  determined, 

the  various  exponents,  n^,  may  be  calculated  by  substituting  a force  or 


moment  value  for  an  angle  in  the  range  which  is  to  be  most  accurately 
approximated  and  solving  for  the  exponent.  The  moment  coefficients  are 
assumed  to  be  linear  over  a +10°  range  about  zero  to  preclude  problems  with 
the  infinite  slope  characteristic  of  (sin  x)n  at  the  origin  for  n < 1. 

Figure  3 presents  one  possible  approximation  to  the  XV-6A  data.  The 
corresponding  functional  expressions  are  included  in  Table  II. 

Using  the  functions  of  Table  I,  the  fuselage  force  and  moment  coefficients 
are  calculated  for  a specific  a and  ay.  They  must  then  be  dimensionalized, 
resolved  into  the  proper  body  axes  and  transferred  to  the  vehicle  CG. 


The  coefficients  are  dimensionalized  using  the  component  of  dynamic 
pressure  in  the  coordinate  plane  corresponding  to  that  of  the  coefficient. 
That  is: 


A = qQ  (A/q0(at)) 
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TABLE  II 

XV-6A  Fuselage  Aero  Approximation 


M/q, 


N/q, 


•C; 

•i; 


A/q  - 0.123  1 

O 1 

COB  « | 1‘8519 

S/q  - 17.89  1 

Bln  « I1-7676 

o 1 

y< 

«i\  - 17-12  | 

Bln  (<v-2)  | 1* 

; 1 Bin  180<*-2> 

0.8046 

1 83 

, I . 180(«-85) 

1 | .1. TJ7  ~ 

' 0.3468 

0.7983 


104 1 ■*»  nsfr . 

099 1 i °-6081 

110 


-8l£«£85 
«<-81  ; or>85 


-70£o-  £70 

y 

» <-70  ; of  >70 

y y 
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S "*  qXy  (b/ (ay)  ) 
Nf  = qxz  (Nf/q0(a)) 


M = qxz  (M/q0(a)> 

N = qxy  (N/q0(ay>)  (3) 


where: 

-1  ,,  2 ^ 2,1/2  . , 
at  = tan  ((v  + w ) / u) 

q^  = 1/2  p (u2  + v2) 

q0  = 1/2  p (u2  + v2  + w2) 

qxz  = 1/2  p (u2  + w2)  (4) 

The  forces  and  moments  acting  in  the  body  axis  directions  at  the  fuselage 
aerodynamic  reference  point  are  given  by: 

Xp  ~ —A  Lp  = 0 

Yp  = S Mp  = M 

ZF  = -Nf  Np  = N (5) 


Finally,  after  transferring  to  the  vehicle  CG: 

AXp  = Xp  ALp  ~ Zp(BLp  ~ BLqq)  — Yp/WL^g  - WLp) 


GYp  = Yf 


AMp  55  Mp  + Xp(WLgg  — WLp)  — Zp(FSgg  — FSp) 


AZp  = Zp 


ANF  = Np  + Yp(FSCG  - FSp)  - Xp(BLF  - BLCG) 

(6) 


Lifting  Surface  Aerodynamics 

Basic  lift  and  drag  characteristics  of  the  three  primary  lifting 
surfaces  (wing,  horizontal  stabilizer  and  vertical  stabilizer)  are  calculated 
as  functions  of  angle  of  attack.  The  approach  used  applies  to  all  three  and 
is  presented  generally.  Particular  effects  which  are  unique  to  only  one 
lifting  surface  are  presented  following  the  general  discussion.  Lift  and 
drag  coefficients.  Cl  and  CD  are  calculated  over  the  entire  angle  of  attack 

range  from  -180  to  +180  degrees  using  primarily  geometric  data  with  some 
basic  aerodynamic  data  estimated  based  on  the  methods  of  references  (d) 
through  (f)  or  similar  sources.  The  calculated  forces  are  assumed  to  act 
at  the  CP  (Center  of  Pressure)  of  the  lifting  surface  in  directions  parallel 
and  perpendicular  to  the  local  velocity  (V^)  which  is  deflected  from  the 


13 


NADC-76312-30 


freestream  velocity  (V)  by  the  induced  angle  of  attack,  a^.  Figure  4 

describes  the  required  lifting  surface  geometric  characteristics.  Appendix  D 
of  reference  (f)  provides  a convenient  method  to  determine  the  CP  location  lor 
a given  planform.  Note  that  this  location  must  be  expressed  in  the  vehicle 
FS,  BL,  and  WL  coordinates  and  for  the  wing  must  be  determined  for  each 
half  of  the  planform.  During  the  development  and  validation  of  this  method 
it  was  found  that  using  the  CP  of  only  the  exposed  portion  of  each  planform 
yielded  the  best  results. 

The  model  for  the  three-dimensional  Cl  is  based  on  the  data  presented 

in  Section  4. 1.3. 3 of  reference  (d).  The  basic  formulation  is  repeated 
here  and  the  reader  may  refer  to  reference  (d)  for  a more  detailed  discussion. 
Where  possible,  empirical  data  curves  from  reference  (d)  have  been  approximated 
with  analytical  functions  for  computational  efficiency. 

For  the  purpose  of  lift  calculations  the  angle  of  attack  range  may  be 
divided  into  two  major  subdivisions.  The  first  or  normal  flow  region  is 
defined  as  that  region  where  the  axial  velocity  component  is  in  the 
negative  X-direction  (-90°  < a < 90°);  the  second  or  reversed  flow  region 
is  defined  as  that  region  where  the  axial  component  of  velocity  is  in  the 
positive  X-direction  (-180°  < a < -90°  and  90°  < a < 180°).  Aerodynamic 
calculations  are  similar  in  both  regions  with  only  the  input  data  being 
different.  Each  of  these  regions  may,  in  turn,  be  divided  into  two 
subregions:  one  for  angles  below  CLmax  (a  < a^)  and  one  for  angles  above 

Climax  (a  > ab).  Here,  angle  of  attack,  a,  is  defined  as  tan  * (w/u)  plus 

the  incidence,  i,  of  the  zero  lift  line  (average  for  twisted  surfaces)  of 
the  surface.  The  following  equations  describe  Cl  for  0°  £ a <_  90°. 

However,  the  equations  hold  for  -90°  <_  a <_  0°  with  a sign  reversal  in  Cl 

and  for  the  reversed  flow  regions  with  the  corresponding  reversed  flow  input 
values  and  the  appropriate  change  in  angle  of  attack  reference. 


In  general,  c£  is  defined  by: 

C^  = CLa(sin  2 a cos  a/2)  + C^^  sin  a cos  alsin  al  (7) 

where: 

CLa  = KW(b)  (CLa)e  (Se/S)  (8) 

KW(b)  « 0.527(1  + d/b)1'^4  + 0.473  (9) 

(CLa)e  ' 2 it  ARe/(2  +\|(2trARe/a0)2(l  + tan2A1/2)  + 4)  (10) 

tan  A1/2  - tan  Ai/a  -(4/AR)  (0. 25 (1-A) / (1  +A))  (11) 


Equation  (9)  represents  the  data  of  Figure  4.3.1.2-10  of  reference  (d). 

The  nonlinear  coefficient,  Cfjaa,  in  equation  (7)  is  calculated  differently 
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for  the  two  angle  of  attack  regions.  For  a _< 

Cv  = (Cm  ) + ACm 

waa  waa  o Naa 

where: 


(CLmax/cOS  “b  ~ CLa  (sin  2 ub)/2) 
(CNaa)o  = (sin  at,  | sin  | ) 


(13) 


2 4 

ACv,  = K cos  ' ( Titan  a/ (2  tan  a,  )) 

Wnim  D 


(14a) 


1/2 

628  J for  J > 0 
22  J for  J < 0 


(1+C1)(1+C2)  - 

(1+C2)AR  tan 

3 

7 



— — J 

— 

J = 0.3(1+C^)AR  cos  A^e 


C1  = 4.470  X3  - 8.125  X2  + 3.712  X - 0.029 


(14b) 

(15) 

(16a) 


C2  = 2.943  X3  - 7.208  X2  + 5.199  X - 0.113  (16b) 

Equations  (14)  and  (16)  are  functional  representations  of  Figures  4.1.3.3-55a 
and  4.1.3.4-24b  of  reference  (d).  For 


CN«a  = (CNaa)o  + ((CNaa)90  " (CNaa)0)(1 


tan  ab  CLa 

7~n~-  ) + D(  -=-!  ) (17) 

tan  a Z . J 


where  (CNaa)0  ls  8tven  by  equation  (13),  (cNaa)90  = !-16  and ; 


D - -1.55  sin  (u(l  - 0.6  ( 


tan  a.  tan  a.  2 

^ ) -0.4  (- )) 

tan  a ' tan  a 


(18) 


Equation  (18)  represents  the  data  of  Figure  4.1.3.3-55a  of  reference  (d). 
The  input  values  of  and  in  the  preceeding  equations  are  determined 

using  the  methods  of  Section  4. 1.3.4  of  reference  (d). 

The  lifting  surface  drag  calculation  is  based  on  the  formulation 
presented  in  reference  (e)  and  is  repeated  here  for  completeness. 

For  a < a. : 


CD  “ CD0  + cDa  a + cDa2  a2 


(19) 
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Remember  that  this  drag  is  parallel  to  the  local  velocity  vector  and  therefore 
does  not  include  induced  drag.  CpQ  is  estimated  using  Section  4. 1.5.1  of 
reference  (d).  Cp^  and  Cpa2  are  normally  zero  for  a < but  may  be  included 

if  available  data  indicates  non-zero  values.  For  a > a,  : 

b 


, (u-TT/2)2(CDb  - 1.2) 

CD  = 

(®b  ~ tr/2  r 


+ 1.2 


(20) 


where : 


cDb  = Cdo  + CDa  ab  + CDa2 


(21) 


F.quations  (7),  (19)  and  (21)  represent  lift  and  drag  referenced  to  the 
direction  of  the  local  velocity  which  must  be  transformed  through  the  induced 
angle  of  attack  to  yield  Cp  and  Cp  defined  in  the  conventional  manner,  i.e., 
perpendicular  and  parallel  to  the  free  stream  velocity  vector  in  the  x-z 
plane.  Performing  this  transformation: 


I t 


CL  = CL  cos  ~ Cp  sin 

(22) 

» » 

Cp  = Cp  sin  + Cp  cos 

(23) 

where : 

= C^/ ("AR  e) 

(24) 

e = 0.527  + 0.1494  AR  - 0.01429  AR2  (25) 

Equation  (25)  is  a functional  description  of  the  data  of  Figure  2.4  of 
reference  (e)  which  is  valid  for  2 < AR  < 7.  Figure  5 presents  Cp  and  Cp 

for  the  XV-6A  Kestrel  wing  calculated  in  this  manner. 

For  a lifting  surface  which  does  not  have  a symmetrical  airfoil  section 
(i.e.,  one  which  produces  a non-zero  pitching  moment  at  the  zero  lift 
condition)  provision  has  been  made  to  include  Cmo  in  the  aerodynamic  model. 

The  value  of  Cmo  for  a particular  geometry  may  be  estimated  using  the  methods 

of  Section  4. 1.4.1  of  reference  (d).  Pitch  moment  due  to  lift  and  drag  at 
angle  of  attack  is  accounted  for  in  the  aero  force  transformation  to  the 
vehicle  CG. 

This  completes  the  general  description  of  the  lift  and  drag  model. 
Additional  effects  for  each  of  the  surfaces  will  now  be  discussed. 

1 . Wing 

For  the  wing,  two  additional  aerodynamic  effects  are  modelled:  the 
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6f  = 0° 


I I I I | 


-180  -90  0 90  180 

a,  degrees 

CD 


FIGURE  5:  XV-6A  CL  AND  CD  APPROXIMATION 
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effect  of  hlgh-lift  flaps;  and  the  effect  of  sideslip,  roll  rate  and  yaw 
rate.  Deflecting  a high-lift  flap  affects  lift  and  drag  at  zero  alpha, 
maximum  lift  and  pitching  moment.  These  affects  are  added  to  the  basic 
aero  model  of  the  wing  as  a function  of  flap  deflection. 

CLr 


CL„ 


" ( 77  )6f 


(26a) 


ACr 


ACDQ 

)6f 


(26b) 


>*f 


(26c) 


tcm  = Cm{ 


(26d) 


Section  6.1  reference  (d)  provides  methods  for  estimating  the  coefficients 
of  equations  (26).  Typical  flap  effects,  as  modelled,  are  shown  on  Figure  5 
for  a flap  deflection  of  60°. 


Increments  in  wing  roll  and  yaw  moment,  1^  and  Nw  (in  stability  axes), 

due  to  sideslip,  3;  roll  rate,  p;  and  yaw  rate,  r;  are  calculated  using 
derivatives  which  are  estimated  by  the  methods  of  reference  (d) , (f),  or 
(g)  and  the  following  equations. 

' 2 r ACle  bw  AClr  ~ 

■v  ■ 1'2pW4  [<Cl6o  + cl) 8 ♦ ^ CL  r + Clp  p)J  (27) 


' 2 ACne  2 bW  A^nr 

N = 1/2  pS  b <z  (C  + —f  CL)B  + 2V“  ( — T 

cf  xz  c : 


2 ACnr 
CL  + -C7-  -D 


Cn)r 


W 


AC. 


2V 


ti£  ACn. 


xz 


‘ ~oC  ci + ^ c”»)p_l 


(28) 


where : 


I 2 2 

vxz 

-kJu  + w 

(29a) 

6 - 

tan  1 (v/Vxz) 

(29b) 

and  Cl  and  Cp  are  the  averages  of  the  left  and  right  wing  values. 
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Now,  dimensionalizing  the  lift,  drag,  and  pitch  moment  coefficients, 
and  transforming  the  forces  and  moments  into  body  axis  directions  yields: 


Hlu.ru  - 1/4»vLsu 

j^L  sin  (a-i)  - CD  cos(a-iTjLWjRW 

(30a; 

ZLU , RU  ■ 1/“>VLSW 

-CL  cos  “ CD  sln  (-iT]LW,RW 

(30b) 

”w  1/2pVxzSWCW  Cmo  + 6fjw 

(30c) 

f 

Ly  = cos  (a-i)w  ■ 

- Nw  sin  (a-i)w 

(30d) 

f 

Nw  = Kf  sin  (a_1)w  ' 

t 

h Nw  cos  (a-i)w 

(30e) 

Transferring  the  forces  and  moments  to  the  vehicle  CG  yields: 

AXW  * *LW  + 

4Lu  ■ Lu  + 2lw(bllu-  blcg)+zru<blru-blcc) 

‘V  ° "V  ' VWmCG-WLlU>-2I.U<FSCG-FSLW> 

",tRW<“-CG-"LRU>-2RW(FSCG-FSRW> 

‘2U  ‘ ZLW  + ZRW  SNW  ’ VXLW<BlLW-BLCG,-XRW(BLRW-BLCG)  <31> 

where  (FS,  BL,  WL)^  and  (FS,  BL,  WL)R^  locate  the  CP's  of  the  left  and 
right  wing  panels. 

2 . Horizontal  Stabilizer 

The  aerodynamics  of  the  horizontal  stabilizer  are  affected  by  the  wing 
dovnwash  field  and  the  vehicle  angular  rates.  Total  effective  angle  of 
attack  is  assumed  to  be: 

<*H  - tan  1 (Wn/UH>  - e + AaR  + iH  (32) 


where: 

“H  ■ “ + ’(“lCG-“lH>  -r  (BLH-BLCG> 

(33a) 

(33b) 
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e " <e/(VCLW 


AaH  - (-/CLw)CLa  aw( 


FVFSw 


w 


(3A) 


(35) 


Equations  (33)  account  for  induced  velocities  at  the  tail  due  to  vehicle 
angular  rates.  Equations  (34)  and  (35)  represent  corrections  for  wing 
downwash  and  lag  of  downwash,  respectively.  For  a canard  configuration, 
these  terms  would  represent  an  upwash  (i.e.,  e/Cj^  < 0).  The  wing  wake 

also  reduces  the  dynamic  pressure  at  the  tail.  This  effect  is  approximated 
by: 

Aqo  2 2 

(q0)H  * (1-  — ) 1/2P  (i^  + wH)  (36) 


where : 


90 


2.42  Cn  V- 
J->o  2 

x' 

£r-  + 0.30 


2 

cos 


(nz  /2zw) 


(37) 


zw  = 0.68  [cDo(  f + 0.15>]  l72 

) 

Equation  (38)  for  2 is  the  semi-thickness  of  the  wake.  The  remaining 
geometry  is  defined  in  Figure  6.  This  formulation  is  based  on  that 
presented  in  Section  4.4  of  reference  (d). 


Body  axis  components  of  the  horizontal  stabilizer  forces  and  moments 
at  its  CP,  located  at  (FS,  BL,  WL)^,  are  given  by: 


*H  ' <"°>H  SH 
ZH  = (qo)H  SH 


|~Cl  sin(a-i)  - Cp  cos  (a-iTj  H 
Pl  cos  (a-i)  - CD  sin  (a-i)]  H 


(39a) 

(39b) 


MH  - (qo)H  SH  ch(CIBo)h  (39c) 

Referenced  to  the  vehicle  CG: 

axh  - *h  alh  - *h(blh-blcg> 

ayh  * 0 AMj,  = mh+xh(wi.cc-wlh)-zh(fscg-fsh)  (40) 
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AZH  = ZH  anh  ' xH^BLH-BLrjG) 

3 . Vertical  Stabilizer 

Angle  of  attack  of  the  vertical  stabilizer  is  defined  by: 

= tan  ^(-Vy/uy)  (41) 

where: 

uv  = u+q(WLCG-WLv)-r(BLv-BLCG)  (42a) 


vy  = v(l-Kv)-p(WI.CG-WLv)+r(FSCG-FSv)  (42b) 

The  coefficient  Ky  in  equation  (42b),  commonly  called  the  sidewash 
coefficient,  accounts  for  the  change  in  lateral  velocity  at  the  vertical 
stabilizer  due  to  flow  around  the  aft  fuselage.  This  coefficient  is 
difficult  to  estimate  and  should  be  set  to  zero  if  applicable  data  are  not 
available.  The  dynamic  pressure  acting  on  the  vertical  stabilizer  is  given 
by: 

(qG)v  = l/2p(uj  + vj)  (43) 

As  for  the  wing  and  horizontal  stabilizer,  the  vertical  stabilizer 
forces  and  moments  are  transformed  to  body  axis  components  and  translated 
to  the  vehicle  CG  resulting  in  equations  (44)  and  (45). 


Xv  = 

<%>VSV 

Cl  sin 

(a-i) 

- Cp  cos 

(a-i) 

V 

(44a) 

Yv  - 

<<Jo)vSv 

Cl  cos 

(a-i) 

+ Cp  sin 

(a-i) 

V 

(44b) 

Nv  - 

^o)vSVCV^Cmo^v 

(44c) 

AXy 

= *v 

ALv  = 

-V'W 

AYV  ■ 

- Yv 

AMy  = 

xv(ULcc' 

-WLV) 

AZv 

= 0 

ANV  = 

"v  + Yv 

(fscg-fsv)-xv<blv-blcg) 

(45) 

Control  surface  effects  have  not  been  discussed  but  are  simulated  by  an 
appropriate  increment  in  the  surface  angle  of  attack.  The  manner  in  which 
this  increment  is  calculated  is  described  in  a later  section  of  this  report. 
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Propulsion-Induced  Aerodyna.n ics 


There  are  at  least  two  significant  propulsion-induced  aerodynamic  etfets 
which  must  be  included  in  the  aerodynamic  model  for  the  vehicle.  The  first 
of  these  is  the  aerodynamic  interference  of  the  propulsion  system  efflux  in 
the  vicinity  of  the  airframe.  This  phenomenon  can  produce  increments  in  the 
vehicle  aerodynamics  which  are  functions  of  engine  thrust,  nozzle  angle  and 
free  stream  velocity  (magnitude  and  direction).  Unfortunately  these  effects 
are  normally  very  configuration  dependent  and  the  available  data  base  is  far 
too  limited  to  permit  complete  modelling  within  the  scope  of  this  effort. 
Review  of  some  available  data  on  several  jet-lift  configurations  has  led, 
however,  to  the  following  general  model  approximations. 

2 3 

= A|  + (B^u  + CjU  + D^u  )sin  0j  (4b) 


(A2u  + B2u2)  (47) 


AL 

T 


(AjV  + 


2 

B3vZ) 


(48) 


AM 

T 


2 3 

A.+(B.u  + C.u  + D.u  )sin  0T 
4 4 4 4 J 


(49) 


The  coefficients  in  equations  (46)  through  (49)  may  be  derived  from  existing 
data  on  the  subject  configuration  or  similar  configurations.  Note  that 
equations  (48)  and  (49)  represent  roll  and  pitch  moment  about  the  vehicle 
CG.  Transforming  to  body  axis  increments  yields: 


AXj  = T ( sin  a-  ^ cos  a)  AL^  = T(  njjr  ) 


AYj  =0  AMj  = T(  —■  ) 


AZj  = T(  — — cos  a — — sin  a)  Nj  = 0 (50) 

Considerable  effort,  both  theoretical  and  semi-empirical,  is  currently  being 
expended  in  this  field  (both  in  and  out  of  ground  effect)  And  this  portion 
of  the  aerodynamic  model  should  be  updated  as  additional  data  become  available. 

The  second  interference  effect  is  that  produced  by  the  engine  (or  fan) 
inlet  momentum  (ram  drag).  For  CTOL  configurations,  inlet  momentum  primarily 
reduces  the  installed  net  thrust  due  to  the  ram  drag  effect.  However,  for 
V/STOL  configurations,  other  significant  forces  and  moments  may  also  be 
produced  due  to  unconventional  inlet  locations  and  large  angles  of  attack 
and  sideslip.  Reference  (h)  indicates  that  a good  approximation  to  inlet- 
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momentum- induced  forces  and  moments  is  achieved  by  applying  the  inlet 
momentum  drag,  w^V,  in  the  direction  of  the  free  stream  velocity  at  a point 

one  inlet  diameter  ahead  of  the  inlet  (along  the  inlet  axis).  With  this 
assumption,  the  following  force  and  moment  equations  were  formulated. 


AX. 


RAM 


w u 
a 


alram  ayram(wlram  «lcg) 


AY 


RAM 


W V 

a 


amram 


-axram(wlram-wlcg>-azram(fscg-fsram^ 


AZ 


RAM 


w w 
a 

g 


anram  ayram(fscg"1'sram) 


(51) 


The  total  weight  flow  of  air  through  the  inlet  is  wa  and  is  approximately 
proportional  to  thrust, 
drag  application  points  for  all  inlets. 


and  WL„,,,  locate  the  centroid  of  ram 
RAM  RAM 


MASS  AND  INERTIAL  CHARACTERISTICS 


Vehicle  gross  weight,  W,  and  moments  of  inertia,  Ix,  Iy,  Iz,  and  Ixz, 

are  required  data  for  the  model.  The  inertias  are  referenced  to  the 
aircraft  body  axis  system  and  the  weight  is  resolved  into  this  system  in 
the  standard  manner. 

AX™  = -W  sin  6 


AY,_  = W cos  0 sin  <p 

W l 


AZ™  = W cos  0 cos  <p  (52) 

0 and  # are  the  aircraft  Euler  angles  in  pitch  and  roll.  The  aircraft  weight 
is  assumed  to  be  constant. 

PROPULSION  SYSTEM  CHARACTERISTICS 

Both  fixed  and  vectorable  nozzle  engines  are  included  in  the  propulsion 
section  of  the  model.  These  nozzles  are  assumed  to  provide  the  primary 
horizontal  and  vertical  thrust  of  the  vehicle.  If,  for  a given  configuration, 
any  of  these  nozzles  are  also  used  for  active  moment  control,  the  moment  control 
is  modelled  by  placing  reaction  control  jets  at  the  same  location.  These  jets 
produce  only  the  incremental  thrust  for  moment  control  and  their  simulation 
is  described  in  a later  section  of  this  report. 

Fixed  Nozzle  Engines 

One  or  two  fixed  nozzles  may  be  simulated.  If  two  are  simulated,  they 
are  assumed  to  be  symmetrically  located  about  the  aircraft  x-z  plane,  as  is 
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normally  the  case.  The  thrust  vector  location  is  described  by  the  FS , 31., 
and  WL  of  Its  application  point;  Its  direction  Is  described  by  yaw,  </,  and 
pitch,  0,  relative  to  the  x-axis  (Figure  7).  Body  axis  components  of  the 
thrust  of  each  jet  are  defined  by  the  standard  axis  transformation  of 
reference  (1). 

AX  = T cos  0 cos  1 (/ 
r J r J 


AYpj  = Tpj  cos  0 sin  i/( 


azfj  “ “tfj  sln  6 


(53) 


Moments  about  the  CG  are  given  by: 


SLfj  - aZrJ(BLFJ-BLcc>-4YFJ(WLCG-WLrj) 


iMFJ  ' 4XFJ(Wl-CC-HLFJ>-aZFJ<rSCG-FSFJ> 


“fj  ' bffj  ^fbcg-fbfj^-bxfj^bbfj-bfcg^  ,W 

Vectorable  Nozzle  Engines 

As  many  as  six  vectorable  nozzles  may  be  simulated  within  the  model  wtiich 
has  been  developed.  Location  is  again  described  by  the  FS,  BL,  and  WL  of 
the  thrust  vector  application  point.  Two  distinct  means  of  describing  the 
orientation  of  the  nozzle  exhaust  direction  are  provided  so  that  the  plane  of 

motion  of  the  nozzle  (the  plane  containing  0 ) may  be  located  appropriately- 

3q 


1.  Roll,  <t> , about  the  x-axis  followed  by  pitch,  © , about  the 

Jo 

displaced  y-axls; 

2.  Yaw,  ii , about  the  Z-axis  followed  by  pitch,  0 , about  the  displaced 

do 

Y-axis . 


Figure  7 depicts  the  two  possible  orientation  schemes. 

Resolution  of  the  thrust  vector  into  body  axis  components  differs  for 
each  orientation  method.  For  the  first  method: 

AX  - -T  sin  (0  + A0  ) 

J J 


AY . » T cos  (0  + A0  ) sin 

J J0  J 
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FIXED  NOZZLE 


VECTORABLE  NOZZLE 


Method  1 


FIGURE  7:  FIXED  AND  VECTORABLE  NOZZLE  ORIENTATION 
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AZ . = -T  cos  (0  + A0  ) cos  $ (5 5u) 

J JQ  J 

For  the  second  method: 

AXt  = -T  sih  (6  + A0  ) cos 

J Jq  j 


AYt  = -T  sin  (0  + A0  .)  sin  ip 

J Jo  J 

AZj  = -T  cos  (0  + A0T)  ( 55b) 

JQ  J 

In  both  methods,  A0^  is  the  change  in  jet  orientation  from  the  reference 

position,  9jo»  due  to  nozzle  vector  control.  Moments  produced  by  these  jets 

about  the  aircraft  CG  are  calculated  in  the  same  manner  as  those  for  the 
fixed  nozzle  jets.  Thrust  management,  both  magnitude  and  direction,  is 
described  in  the  Control  System  section  of  this  report. 

Angular  Momentum  Effects 

Spinning  masses  within  the  engine  (compressors,  turbines,  etc.)  and 
their  associated  angular  momenta,  H,  couple  with  the  angular  rate  of  motion 
of  the  aircraft  to  produce  increments  in  the  moments  acting  about  the 
\ aircraft  CG.  In  the  case  of  conventional  aircraft,  where  aerodynamic 

damping  moments  are  significant,  these  gyroscipic  moments  are  usually 
neglected.  However,  for  V/STOL  at  low  speed,  these  moments  may  be  quite 
significant.  For  this  reason , these  moments  are  included  in  the  model. 

The  angular  momentum  of  the  fixed  nozzle  engines  is  assumed  to  be  about 
the  axis  defined  by  the  thrust  vector,  with  positive  direction  being 
clockwise  looking  in  the  direction  of  the  thrust.  The  body  axis  components 
are  therefore  defined  by: 


= H cos  0 cos 
Ry  = H cos  0 sin  \p 

Hz  = H sin  0 (56) 

The  moment  produced  by  engine  angular  momentum  is  derived  in  reference  (f)  as: 

ALg  - rRy  - qHz 


amg  = p«z  - rHx 

ANn  = qHx  - PHY  (57) 
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The  angular  momentum  of  the  fixed  nozzle  engines  is  assumed  to  be  proportional 
to  the  thrust  of  the  engine.  If  this  assumption  proves  to  be  too  gross  for 
the  engine  under  consideration,  the  model  described  for  the  vectorable 
nozzle  engines  may  be  used. 

Since  the  angular  momentum  vector  for  a vectorable  nozzle  engine  cannot 
be  assumed  to  coincide  with  the  thrust  vector,  a separate  set  of  orientation 
angles  must  be  defined.  The  angular  momentum  vector  is  located  relative  to 
the  aircraft  x-axis  by  a yaw  and  pitch  angle  in  the  same  manner  as  the  fixed 
nozzle  engine  thrust  vector.  Increments  in  aircraft  moments  are  calculated 
in  the  same  manner  as  that  described  for  the  fixed  nozzle  engines.  The 
magnitude  of  the  angular  momentum  is  assumed  to  be  a function  of  thrust 
defined  by: 


H ( 
max 


A + BT  + CT 
100 


(38) 


2 

The  coefficients  A,  B,  and  C are  in  units  of  %,  %/newton  and  X/newton", 
respectively,  and  must  be  determined  from  available  engine  data.  Angular 
momentum  contributions  are  associated  with  each  nozzle  and  its  respective 
thrust  when  there  is  more  than  one  nozzle  per  engine. 


CONTROL  SYSTEM  CHARACTERISTICS 

The  control  system  model  consists  of  longitudinal  stick,  lateral  stick, 
and  rudder  pedals  for  regulation  of  both  aerodynamic  control  surfaces  and 
reaction  control  nozzles.  In  addition  to  these  conventional  controls,  a 
thrust  management  system  consisting  of  throttles  and  nozzle  controls  is  also 
included  in  the  model.  In  general,  the  control  system  model  relates  a total 
of  19  outputs  to  8 inputs  (Figure  8).  To  facilitate  detailed  description 
of  the  system,  it  is  divided  into  three  segments:  Aerodynamic  control, 

reaction  control  and  thrust  management. 

Aerodynamic  Control  Surfaces 

The  effect  of  deflecting  an  aerodynamic  control  surface  is  modelled  as 

an  incremental  change  in  angle  of  attack  of  the  lifting  surface  of  which  the 

control  surface  is  an  integral  part.  For  example,  rudder  deflection  is 

simulated  as  an  incremental  change  in  vertical  stabilizer  angle  of  attack, 

Aay.  The  longitudinal  stick  regulates  horizontal  stabilizer  angle  of  attack, 

the  lateral  stick  regulates  left  and  right  wing  angle  of  attack  and  the 

rudder  pedals  regulate  vertical  stabilizer  angle  of  attack.  Effective  angle 

of  attack  change  due  to  control  surface  deflection  may  be  estimated  by 

applying  the  methods  of  Section  6 of  reference  (d).  Figure  9 presents  an 

example  of  the  assumed  form  for  control  surface  gearing  and  sign  conventions. 

Note  that  in  the  case  of  lateral  control,  both  Aa  and  Aa.  . are  calculated 

RW  LW 

and  are  assumed  to  be  equal  in  magnitude  but  opposite  in  direction. 

Reaction  Control  Nozzles 


Up  to  a maximum  of  ten  reaction  control  jet  nozzles  are  included  in  the 
model.  Each  nozzle  is  located  in  the  aircraft  axis  system  in  a manner 
identical  to  that  used  for  the  fixed  nozzle  engines  (FS,  BL,  WL,  \p  and  9). 
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o ’ ' 50  ' 100 


Percent  Control 

FIGURE  9:  CONTROL  SURFACE  GEARING  AND  SIGN  CONVENTIONS 
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Forces  and  moments  due  to  the  thrust  produced  by  each  nozzle  are  calculitc! 
using  the  following  expressions. 

fiXRJ  - Tr  cos  0 cos  t a lrj  = azrj^BIjrj-BLCG^~AYRJ^WLCG-WLRJ^ 

4YRJ  ■ TR  COS  9 8l”  * 4MRJ  ■ 


4zrj  ■ -tr  s1"  9 4NR.I  • 4YRJ<FSCC-FSRJ>-SX,U<BLm-BLCG)  <5,) 

Variations  in  RCS  nozzle  thrust,  T^,  due  to  control  and  main  engine  bleed-air 
availability  are  described  in  the  following. 

As  described  above,  each  nozzle  is  oriented  according  to  its  positive 
thrust  direction;  however,  for  control  each  nozzle  is  assumed  to  be  capable 
of  thrusting  in  both  positive  and  negative  directions.  The  assumed  variation 
of  nozzle  thrust  with  control  motion  is  shown  in  Figure  10.  The  center 
control  position  with  respect  to  RCS  nozzle  thrust  is  defined  as  <50  with 

a deadband  of  width  6^  on  either  side.  Once  the  control  is  moved  out  of 

the  deadband  the  nozzle  thrust  increases  in  direct  proportion  to  control 
deflection  up  to  maximum  nozzle  thrust,  Tj^  . This  increase  takes  place 

over  the  range  of  control  motion,  The  maximum  thrust  may  be  different 

for  positive  and  negative  control  deflections.  A nozzle  which  thrusts  in 
only  one  direction  is  simulated  by  setting  the  maximum  thrust  in  the  opposite 
direction  to  zero.  Each  nozzle  thrust  is  calculated  in  this  manner  and  then 
all  are  summed  to  obtain  total  RCS  forces  and  moments. 

There  may  be  a dynamic  lag  in  the  response  of  the  RCS  nozzle  thrust  to 
control  command.  The  possibility  of  up  to  a second  order  lag  is  accounted 
for  within  the  RCS  model.  In  La  Place  Transform  notation,  the  expression 
for  a second  order  lag  is: 


T 1 

Tc  = ( r x s+ 1 ) (t2s+1)  (6° 

Tc  is  the  commanded  thrust  as  determined  by  the  method  of  the  preceeding 
paragraph  and  and  t ^ are  time  constants  in  seconds.  A first  order  lag  is 

modelled  by  setting  X2  to  zero  in  equation  (60).  The  differential  equation 

form  of  equation  (60)  is: 


txT  + T - Tc 


T1T2T  + (ti+t2)T  + T * Tc 


(First  Order) 
(Second  Order) 


(61a) 

(61b) 
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Normally  the  RCS  nozzles  are  powered  by  air  which  is  bled  from  one  r 
mere  of  the  primary  engines.  If  this  is  the  case,  the  thrust  of  the  RCS 
nozzles  may  vary  with  engire  thrust.  This  effect  is  approximated  by: 

* t t A(ET/1000)  +B  (ET/1000)2  .... 

(tr)  , tr(  Too 

K actual 

where  ET  is  the  sura  of  the  thrusts  of  all  engines  which  affect  the  RCS  and 
A and  B are  input  coefficients  in  units  of  %/l(P  newtons  and  %/10^  newtons, 
respectively. 

Thrust  Management 

The  thrust  management  system  consists  of  throttle  control  on  both  fixed 
and  vectorable  nozzle  engines  and  thrust  vector  (nozzle  angle)  control  on 
the  vectorable  nozzle  engines.  For  the  purposes  of  thrust  management,  the 
vectorable  nozzle  engines  are  divided  into  two  groups,  each  with  its  own 
throttle  and  vector  control.  This  grouping  facilitates  the  simulation  of 
interconnects  between  throttles  or  nozzle  controls  of  one  set  of  nozzles 
and  the  other. 

Fixed  nozzle  thrust  is  assumed  to  vary  linearly  with  respect  to  the 
throttle. 


tfj  *tfW5t 


(63a) 


Thrust  of  the  vectorable  nozzle  engines  is  also  assumed  to  vary  linearly 
with  throttle  deflection. 


(tV'Vi6t, 


rj2  = (VV26T2 


(63b) 


For  the  sake  of  generality,  nozzle  angle  is  assumed  to  he  a function  of 

both  throttle,  6_,  and  vector  control,  S„. 

I 0 


■ <svvi\ + vv 

“J2  • <*W2st2  * f2<Je2> 


(64a) 


(64b) 


Functions  f^  and  f2  are  assumed  to  be  piecewise  linear  functions  of 
and  6q2>  respectively. 

To  allow  for  interaction  between  vectorable  nozzle  control  sets  1 and 
2,  the  following  functional  relations  are  assumed. 

= A_  + B 6 + C (65a) 

2 1 1 
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se2  = A9  + 


Bel501  + ce'se] 


(65b) 


Setting  all  coefficients  to  zero  in  either  of  equations  (65)  indicates  that 
there  is  no  interaction  between  the  control  sets.  A typical  nozzle  angle 
interconnect  schedule  and  the  means  of  modelling  it  are  shown  in  Figure  11. 


ANALYSIS 

Once  the  incremental  forces  and  moments  being  generated  by  each  of  the 
model  components  are  calculated  as  functions  of  the  total  vehicle  state, 
they  are  appropriately  summed  and  used  in  the  following  analyses.  The 
analyses  described  consist  of: 

1.  Non-linear  trim  calculation; 

2.  Stability/control  derivative  estimation; 

3.  Small  perturbation  stability  analysis; 

4.  Non-linear  time  response  calculation; 

5.  Least-squares  vector  analysis. 

NON-LINEAR  TRIM  CALCULATION 


A Newton-Raphson  iteration  scheme  is  used  to  determine  the  required 
control  inputs  and  vehicle  attitudes  for  a preselected  trim  state.  Since 
the  remainder  of  the  analysis  requires  the  vehicle  to  initially  be  in  a 
trim  condition,  this  calculation  is  essential.  The  iteration  procedure 
is  described  first  followed  by  definition  of  various  possible  categories 
of  trim. 


Fi  = fi(xj  * V : 


The  force  and  moment  equations  for  the  vehicle  may  be  expressed  by: 

i = 1,  6 
j = 1,  11 

k = 1,  6 (66) 

where  the  are  the  total  forces  and  moments  calculated  using  the  previously 

described  model,  x^  is  a vector  consisting  of  control  deflections  and 

vehicle  attitudes  and  is  a vector  of  vehicle  linear  and  angular  velocities. 
The  vector  x.  consists  of  the  following  elements. 

1.  Fixed  nozzle  engine  throttle,  6 


2.  Longitudinal  stick,  <$ 
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3. 

Lateral  stick,  6 

A. 

Rudder  pedals,  SD 
K 

5. 

Yaw  attitude, 

6. 

Pitch  attitude,  0 

7. 

Roll  attitude,  <p 

8. 

Vectorable  nozzle  engine 

throttle  1,  <S 

T1 

9. 

Vectorable  nozzle  engine 

throttle  2,  6_ 

2 

10. 

Nozzle  vector  control  1, 

\ 

11. 

Nozzle  vector  control  2, 

s 

The 

)T, 

desired  trim  condition  is 
five  of  the  eleven  control 

specified  by  trim  forces  and  moments, 
variables,  x^ , and  the  six  velocity 

components,  V^.  Care  must  be  exercised  when  selecting  the  six  control 

variables  to  be  used  for  trim  so  that  control  is  provided  for  each  of  the 
six  degrees  of  freedom.  With  the  above  specification,  the  system  to  be 
solved  for  trim  reduces  to: 


(Fi)T  = fi^x^  5 


i = 1 , 6 
j - 1,  6 


(67) 


The  trim  solution  begins  with  an  initial  estimate  of  the  control  vector, 
x . To  satisfy  the  specified  trim  condition: 


(F  ) - f (x  + Ax  ) = 0 

11  1 1 o 1 


(68) 


Therefore  the  trim  problem  reduces  to  solving  equation  (68)  for  Ax..  This 


is  accomplished  by  an  iterative  process  in  the  following  manner. 

fi(x  + Ax.)  is  expanded  with  second  and  higher  order  terms  neglected 
Jo  J 
yield ing : 


(Fi)T 


fj(x,  ) 
1 jo 


ifJ 

Ox , 


)Ax.  = 0 
x,  j 
jo 


(69) 


Solving  for  Ax 


j‘ 
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Ax 


j 


((Fi)T 

V 


f . (x.  )) 
*•  JO  y 


AF 


i 


(70) 


a 


difference 


is  a six-by-six  matrix 
technique 


evaluated 


at 


x.  and 
jo 


is  calculated  using 


f . (x . + 6x . ) - f . (x . ) 

1 JQ  J 1 JQ 
6x . 

J 


(71) 


The  magnitude  of  6x.  is  assumed  to  be  1/2  of  the  maximum  allowable  change 

in  Ax.  per  iteration.  Note  that  if  the  matrix  of  equation  (71)  is  singular 

(inverse  does  not  exist)  a trim  solution  cannot  be  obtained  using  the  six 
selected  control  variables. 

After  Ax.  is  calculated,  a new  control  vector  estimate,  x.  = x.  + Ax. 


J 


jo 


J 


is  made.  The  increment  Ax.  is  limited  by  a preselected  maximum  value  which 

decreases  as  AF^  approaches  zero  (as  a trim  solution  is  approached).  This 

improves  the  convergence  of  the  iteration  scheme  when  severe  nonlinearities 
are  present.  The  process  is  repeated  until  AF^  is  within  allowable  limits. 

At  this  point,  the  current  value  of  Xj  is  the  trim  solution  for  the  control 

vector. 


The  allowable  trim  conditions  which  have  been  considered  may  be  grouped 
into  three  categories: 

1.  Trim  with  zero  angular  rate; 

2.  Push-overs  and  pull-ups; 

3.  Coordinated  turns. 

The  required  trim  forces  and  moments,  (F^)  , are  calculated  by  substituting 
the  appropriate  velocities  and  accelerations  into  equations  (72). 

X = mu  + qw-rv  (72a) 


Y 


T 


mv  + ru-pw 


(72b) 
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ZT 

mw  + pv-qu 

(72c) 

L_  = 

qr ( I -I  ) - 

pq  i 

(72d) 

T 

z y 

xz 

*T  " 

Pr«x-i2)  + 

2 2 

(P  -r  >IXZ 

(72e) 

nt  = 

"W  - 

^ Txz 

(72f) 

STABILITY /CONTROL  DERIVATIVE  ESTIMATION 

After  trim  has  been  established,  the  aircraft  stability  derivatives  may 
be  evaluated  at  the  trim  point.  These  derivatives  indicate  the  sensitivity 
of  the  aircraft  body  axis  forces  and  moments  to  incremental  changes  in  the 
body  axis  linear  and  angular  velocities.  Each  velocity  is  incremented  and 
the  resulting  changes  in  the  forces  and  moments  are  calculated.  The  stability 
derivatives  are  formed  by  dividing  the  force  and  moment  increments  by  the 
corresponding  velocity  increments. 

The  accuracy  of  this  calculation  is  directly  dependent  upon  the  value 
chosen  for  the  velocity  increments.  Too  small  a value  might  result  in 
errors  due  to  small  differences  in  large  numbers;  too  large  a value  might 
overlook  the  effect  of  localized  non-linearities.  A good  rule  of  thumb  to 
use  in  determining  the  increment  is  to  use  2 percent  of  the  trim  velocity 
for  linear  velocity  increments  and  0.3  percent  of  the  trim  velocity  (in 
metres/second)  for  angular  velocity  increments  (in  units  of  radians  per  second). 
These  increments  result  in  approximately  1 degree  increments  in  the 
individual  surface  angles  of  attack. 

The  control  derivatives  are  calculated  using  equation  (71)  in  the  same 
manner  used  during  the  trim  iterations. 

SMALL  PERTURBATION  STABILITY  ANALYSIS 

The  derivatives  described  in  the  preceeding  paragraphs  may  be  used  in  a 

linearized,  small  perturbation  aircraft  stability  analysis.  The  present 

method  uses  the  standard  3 DOF  uncoupled  analysis  as  developed  in  Flight 

Dynamics  references  such  as  reference  (f).  The  equations  which  form  the 

basis  of  the  analysis  are  presented  in  Table  III.  These  equations  are 

formulated  in  a stability  axis  reference  frame  with  the  controls  fixed. 

Trim  sideslip  and  flight  path  angle  are  assumed  to  be  zero  and,  therefore, 

these  equations  are  valid  only  when  these  assumptions  are  satisfied  by  the 

trim  condition.  The  acceleration  derivatives  (Z.  and  M.)  appearing  in 

w w 

Table  III  are  defined  by: 

Zi"  1/2>VFVFSCC)(C^H<Cl«)W<e/<V  <”*> 


t 

M. 

w 


(FVFSCC)Zw 


(73b) 
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The  calculated  stability  derivatives  must  be  transformed  to  stability  axes 
for  use  in  the  equations  of  Table  III.  This  is  easily  accomplished  using 
the  method  of  Appendix  B of  reference  (f)  or  another  similar  reference. 

The  characteristics  equations  for  both  longitudinal  and  lateral  degrees 
of  freedom  are  formulated  from  the  equations  of  Table  III  and  solved  for 
the  characteristic  roots.  Common  control  transfer  functions  may  also  be 
calculated.  This  procedure  has  been  automated  and  is  described  in  Volume  II 
of  this  report. 

NON-LINEAR  TIME  RESPONSE  CALCULATION 

Once  trim  has  been  established,  the  aircraft  model  may  be  "flown" 
through  a time  dependent  trajectory  which  is  driven  by  a combination  of  input 
forcing  functions  (e.g.,  control  motion,  gust  disturbance,  etc.).  The 
maneuver  time  history  is  computed  by  integrating  the  fully  coupled,  six 
DOE  equations  of  motion  of  the  vehicle  (Table  IV);  recomputing  the  vehicle 
non-linear  forces  and  moments  at  each  time  step.  A four-cycle  Runge-Kutta 
technique  is  used  to  perform  the  necessary  integrations. 

LEAST-SQUARES  VECTOR  ANALYSIS 

A simplified,  least-squares  trigonometric  curve  fit  is  included  in  the 
programmed  analyses  to  facilitate  amplitude  and  phase  angle  analyses  of 
oscillatory  vehicle  behavior.  Portions  of  the  time  history  output  of  the 
program  are  fit  using  general  equations  of  the  form: 

x = C^  sin  ait  + C 2 cos  cot  + C^  (74) 

Coefficients  C^,  C^,  and  C^  are  determined  by  solving  the  following  three 
simultaneous  equations  when  N data  points  of  the  form  (x.,  t^)  are  available 
and  frequency,  to,  is  assumed. 

£xi  = ClZ  sin  oit^  + C££  cos  cot^  + NC.^  (75a) 

2 

ZXj  sin  ui  tj  = C^L  sin  ut^  + C2I  sin  cot^  cos  ut.  + C^I  sin  ut , (75b) 


2 

Lx.  cos  ut.  = C,£  sin  cot,  coscot,  + C„I  cos  10  t . + C_L  cos  co  t , 

i il  i i2  t3  i 

Amplitude  and  phase  are  determined  by: 

Amplitude  = (C^  + C^)  ^ 

Phase  = tan  ^(C„/Cj) 


(75c) 


(76a) 

(76b) 
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TABLE  IV 

SIX  DEGREE  OF  FREEDOM  EQUATIONS  OF  MOTION 

X 

u ■ — + rv  - qw  - g sin  0 

IB 

Y 

v-i+pw-ru+g  cos  0 sin  $ 

■ 

Z 

W « — + qu  - pv  + R COS  0 COS  4 

IB 

iz  r iv_ 

P - 2 L + (N  - qr  !xr)  - qr  (I*  - Iy) 

+ pq  Ixr  (1  - Iy-  IgS 

q - |^M  - Ixz  (P2  - r2)  -rp  <IX  - IZ)J 

H 

' Hr*’] 


- I„ 


(L  + pq  I ) - pq  (Iy  - Ix) 


-qr  Ixz  (1 


t ■ (q  sin  <p  + r cos  4>)  sec  0 

0 <■  q cos  0 - r sin  4> 

♦ » p + q sin  $ tan  0 + r cos  * tan  0 


cos  0 cos  \)r 

sin  $ sin  0 cos  * 
-cos  $ sin  ^ 

cos  t sin  0 cos  ^ 
_+sin  $ sin  ^ 


cos  9 sin  i|i 
sin  $ sin  0 sin  $ 

+COS  $ COS  \|| 

cos  9 sin  0 sin  ^ 
-sin  9 cos  iji 


-sin  0 

u 

sin  $ cos  6 

V 

w 

cos  9 cos  0 

L 
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VALIDATION 

In  order  to  validate  the  models  and  methods  described  in  this  report, 
the  XV-6A  Kestrel  (Figure  12)  was  modelled  using  basic  geometry  and 
aerodynamic  data  and  analyzed  using  the  computerized  prediction  method. 

The  results  were  compared  with  existing  wind  tunnel  and  flight  test  data  from 
references  (j),  (k) , and  (1).  In  general,  the  results  correlated  reasonably 
well  with  the  reference  data. 

Appendix  A contains  a description  of  the  model  input  data  used  and  some 
of  the  data  comparisons  resulting  from  the  validation  analysis. 


CONCLUSIONS  AND  RECOMMENDATIONS 

A six  DOF  V/STOL  Stability  and  Control  modelling  and  prediction  method 
has  been  developed,  computerized  and  validated  with  existing  vehicle  wind 
tunnel  and  flight  test  data.  The  method  is  capable  of  analyzing  both 
static  and  dynamic  characteristics  of  jet  lift  and  lift/cruise  fan  V/STOL 
configurations  in  hover  and  transition  flight.  The  method  is  most  suited  to 
a preliminary  design  environment  where  wind  tunnel  data  are  not  available  or 
are  incomplete. 

As  additional  data  become  available  on  V/STOL  configurations,  further 
validation  should  be  attempted  and  modifications  should  be  made  to  the 
method,  if  warranted,  to  enhance  its  versatility  and  credibility. 
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FIGURE  12:  XV-6A  KESTREL,  3-VIEW 
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LIST  OF  SYMBOLS 


Symbol 

Description 

Units 

ao 

2-D  Lift  Curve  Slope 

1/rad 

A 

Axial  Force 

N (lb) 

B,  C,  D 

Miscellaneous  coefficients 

- 

AR 

2 

Aspect  Ratio,  b /S 

- 

b 

Span 

m (ft) 

BL 

Butt  Line,  positive  on  right  side  of 
aircraft 

cm  (in) 

c 

Mean  Aerodynamic  Chord 

m (ft) 

cr 

Root  Chord 

m (ft) 

ct 

Tip  Chord 

m (ft) 

cD 

Drag  coefficient,  D/q0S,  parallel  to 
free  stream  velocity 

- 

CD 

Drag  coefficient,  D'/q0S,  parallel  to 
local  velocity 

- 

CDo 

Drag  Coefficient  at  Zero  Lift 

- 

% 

Change  in  profile  drag  with  angle  of 
attack 

1/rad 

CDa2 

Change  in  profile  drag  with  angle  of 
attack  squared 

1/rad^ 

CG 

Center  of  Gravity 

- 

% 

Roll  moment  coefficient  due  to  roll 
rate 

1/rad 

Clr 

Roll  moment  coefficient  due  to  yaw 
rate 

1/rad 

ClB 

Roll  moment  coefficient  due  to 
sideslip 

1/rad 

cL 

Lift  coefficient  perpendicular  to  free 
stream  velocity,  Lf/q0S 

- 

» 

cL 

Lift  coefficient  perpendicular  to  local 
velocity,  Lj./q0S 

- 
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Symbol 

Ci 


l-n 


r 


CP 

d/b 

0 


t* 


FS 

g 

H 

i 


Kv 

kW(B) 


LIST  OF  SYMBOLS  ( Continued) 
Descript  ion 

Change  in  lift  coefficient  with 
angle  of  attack 

Zero  lift  pitch  moment  coefficient, 
M/q0ST 

Fitch  moment  due  to  flap  deflection 

Yaw  moment  coefficient  due  to  roll 
rate 

Yaw  moment  coefficient  due  to  yaw 
rate 

Y.iw  moment  coefficient  due  to 
8 Ides  1 ip 

Non-ilnear  normal  force  coefficient 


Units 

1/rad 


1 /deg 
1/rad 

1/rad 

1/rad 

1 / rad“ 


Center  of  aerodynamic  pressure 
Ratio  of  body  diameter  to  span 

Drag  force  parallel  to  free  stream  N (lb) 

velocity,  positive  in  direction  of 

velocity 

Oswald  efficiency  factor  - 

Fuselage  Station,  positive  aft  cm  (in) 

General  force  or  moment  N (lb)  or  N.m  (ft. lb) 


Acceleration  of  gravity 

Angular  momentum 

Lifting  surface  incidence  to 
fuselage  x-axis 

Moments  of  inertia  about  body 
x,  y,  and  z axes 

Product  of  inertia 


kg.m 


m/sec 

2, 

/ sec 


2 


(ft/sec2) 

2 

(slug. ft  /sec) 
deg 


kg.m2  (slug,  ft*") 
kg.m2  (slug. ft2) 


Sidewash  coefficient 


Wing/body  interference  factor 


kl 
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Symbol 
L,  M,  N 

Lf 

m 

MAC 

ni 

Nf 

P,  q,  r 

qo 

s 

S 

t 

T 

T 

*c 

H 

U,  V,  w 

V 

wa 


LIST  OF  SYMBOLS  (Continued) 
Description 

Body  axis  roll,  pitch  and  yaw 
moments 

Lift  force  perpendicular  to  free 
stream  velocity  vector 

Mass 

Mean  Aerodynamic  Chord 

Exponents  in  fuselage  aero 
approximations 

Normal  force,  positive  in  negative 
z-directicn 

Body  axis  roll,  pitch  and  yaw  rates 

2 

Dynamic  pressure,  l/2pV 

La  Place  operator 

Planform  area 

Time 

Thrust 

Commanded  thrust 

Euler  transformation  matrix 
Body  axis  x,  y,  and  z velocities 
Total  velocity 
Inlet  air  weight  flow 


Units 

N .m  ( f t . lb) 

N (lb) 

kg  (slug) 
m (ft) 

N (lb) 

rad/ sec 
N/m2  (lb/ft2) 
1/sec 
m2  (ft2) 
sec 
N (lb) 

N (lb) 

m/sec  (ft/sec) 
m/sec  (ft/sec) 
N/sec  (lb/sec) 


W Weight 

WL  Water  Line,  positive  up 

Xj  Control  vector 

x'  Distance  from  wing  trailing  edge 

to  horizontal  stabilizer  quarter- 
chord,  measured  along  wing  wake 
centerline 


N (lb) 
cm  (in) 


m (ft) 
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LIST  OF  SYMBOLS  (Continued) 


Symbol 

Description 

Units 

X,  Y,  Z 

Forces  in  the  x,  y,  and  z body 
axis  directions 

N (lb) 

zw 

Wing  wake  semi-width  at  horizontal 
stabilizer  location 

m (ft) 

» 

z 

Normal  distance  from  wing  wake 
centerline  to  horizontal  stabilizer 

m (ft) 

l 

Angle  of  attack,  tan  * (w/u) 

deg 

ab 

Angle  of  attack  corresponding  to 
CL 

max 

deg 

ai 

Induced  angle  of  attack 

deg 

aT 

Total  angle  of  attack,  tan  ^(v^+w^)  ^/u) 

deg 

u 

y 

Yaw  angle  of  attack,  tan  ^(v/u) 

deg 

ai*  u2 

Angles  of  attack  corresponding  to 
maximum  fuselage  moments 

deg 

e 

Sideslip  angle,  sin  * (v/V) 

deg 

*'a 

Aileron  deflection,  positive  for  left 
wing  trailing  edge  down 

deg 

4d 

Semi-width  of  RCS  controller  dead 
band 

deg 

5e 

Elevator  deflection,  positive  for 
trailing  edge  down 

deg 

6f 

Wing  flap  deflection 

deg 

«r 

Rudder  deflection,  positive  for 
trailing  edge  left 

deg 

4R 

Rudder  pedal  deflection,  positive  for 
left  pedal 

cm  (in) 

4 RAMI' 

RCS  versus  6 ramp  length 

cm  (in) 

6s 

Longitudinal  stick  deflection,  positive 
for  forward  stick 

cm  (in) 
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Symbol 


P 

Ti 


LIST  OF  SYMBOLS  (Continued) 

Description  Units 

Fixed  nozzle  engine  throttle  deflection  cm  (in) 

Vectorable  nozzle  engine  throttle  cm  (in) 

deflection 

Lateral  stick  deflection,  positive  cm  (in) 

for  right  stick 

RCS  controller  center  position  cm  (in) 

Vectorable  nozzle  angle  controller  cm  (in) 

deflection 

Downwash  angle  deg 

Nozzle  angle  deg 


Taper  ratio,  ct/c 

Sweep  angle  of  the  n*-*1  fractional  chord  deg 

line 

3 3 

Atmospheric  density  kg/m  (slug/ft  ) 

Time  constant  sec 


>>,  $ 

u> 

Subscripts 

CG 

e 

F 

FJ 

G 

H 

I 

J 


Euler  yaw,  pitch,  and  roll  angles  deg 

Frequency  rad/sec 

Referenced  to  vehicle  center  of  gravity 
Property  of  exposed  portion  of  planform 
Fuselage  parameter 
Fixed  nozzle  engine  parameter 
Gyroscopic  (angular  momentum)  effect 
Horizontal  stabilizer  parameter 
Propulsion  induced  aero  effect 
Vectorable  nozzle  engine  parameter 
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Subscripts 

le 

LW 

max 

o 

RAM 
R,  RJ 
RW 
T 

V 
WT 

W 

X 

xy 

xz 

Y 

Z 


LIST  OF  SYMBOLS  (Continued) 
Descript  ion 

Referenced  to  planform  leading  edge 

Left  wing  parameter 

Local  maximum  value 

Reference  condition 

Inlet  momentum  effect 

Reaction  jet  parameter 

Right  wing  parameter 

Trim  condition 

Vertical  stabilizer  parameter 
Weight  component 
Total  wing  parameter 
X body  axis  component 
Component  in  x-y  plane 
Component  in  x-z  plane 
Y body  axis  component 
Z body  axis  component 


Mathemat ical 

A ( ) 

t 

( ) 

( ’),  ( " ) 

Xy 


Symbols 

Incremental  quantity 

Stability  axis  quantity 

First  and  second  time  derivatives 

Partial  derivative  of  X with  respect  to  y 


( 

( 

( 


-1 


) 


Matrix  inverse 
Matrix  transpose 
Quantity  evaluated  at  x 
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To  validate  both  the  modelling  technique  and  the  calculative  procedures 
which  are  described  in  this  report,  an  existing  V/STOL  configuration  was 
modelled  and  the  resulting  predictions  were  compared  with  flight  test  results 
Because  of  data  availability,  the  XV-6A  Kestrel  (Figure  A-l)  was  selected  to 
be  the  sample  configuration.  The  scope  of  the  validation  was  limited  by 
the  available  flight  data.  Physical  characteristics  of  the  XV-6A  were  taken 
from  reference  (j)  and  basic  aerodynamic  data  were  taken  from  reference  (k) . 

The  vehicle  model  input  is  presented  in  five  sections:  (1)  weight  and 

inertia  characteristics;  (2)  fuselage  aerodynamics;  (3)  wing  and  empennage 
aerodynamics;  (4)  propulsion  system  characteristics  and  propulsion  induced 
aerodynamics;  (5)  control  system  characteristics.  Data  for  each  model 
subsection  are  described  in  detail  in  the  following. 

The  nominal  vehicle  weight,  CG,  and  inertias  are  listed  in  Table  A-l. 


TABLE  A- I 

XV-6A  WEIGHT  AND  INERTIAS 
Gross  Weight:  53,379  N 

2 

Inertias:  Ix  = 5016  kg.m 

Iy  = 30,370  kg.m2 
Iz  = 32,946  kg.m2 
Ixz  = 2300  kg.m2 


CG  location  (10.2%tr): 

FS 

556.8 

BL 

0 

WL 

251.0 

Fuselage  aerodynamics  were  estimated  using  the  method  of  reference  (c) 
for  the  complete  angle  of  attack  and  sideslip  range.  The  estimated  aero- 
dynamics were  approximated  with  the  functions  described  in  this  report. 

Figure  A-2  graphically  presents  both  the  estimated  and  curve  fit  data  and 
Table  A-II  lists  the  functional  expressions  which  were  used  for  the  validation 
calculat ions . 

Wing  and  empennage  geometrical  data  and  basic  aerodynamic  characteristics 
are  listed  In  Table  A-III.  The  downwash  and  flap  effects  were  calculated 
using  Sections  4.4  and  6.1  of  reference  (d) . The  lateral/directional  effects 
were  estimated  using  the  charts  and  methods  of  Appendix  B of  reference  (f). 

The  XV-6A  aircraft  is  powered  by  the  Rolls  Royce/Br istol  Pegasus  5 fan 
lift/thrust  engine.  Exhaust  gas  (both  fan  and  turbine)  is  forced  through  four 
swivelling  nozzles — two  on  each  side  of  the  fuselage  producing  approximately 


t,  degrees 


FIGURE  A-2:  XV-6A  FUSELAGE  AERODYNAMIC  DATA 
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TABLE  A- II 

XV- 6 A FUSELAGE  AERO  DATA  FUNCTIONS 


A/qo  ■ 0.123  | cos  a 


S/q  “ 17.89  sin  a 
o Y 


1.8519 

1.7676 


f 9.585  | 

sin 

t-6.533  | 

sin 

f 8.104  | 

sin 

1 -9.099 

| sin 

Nf/qo  - 17.12  | sin  (»-2)  | 

180 (»-2)  | °*8046 

180(«-85)  0,3468 

97 


1.7533 


1800- 


HO 


0.7983 


0.6081 


-81^a-<85 
o-<-81  ; o->85 


-70^»  <70 

y 

o-  <-70  ; o-  >70 

y y 


NADC-76312-30 


WING  AND 

TABLE  A- I II 
EMPENNAGE  AERO 

DATA 

Wing 

Horizontal 

Stabilizer 

Vert  leal 
Stabil izer 

2 

Reference  Area,  S(m  ) 

17.32 

4.41 

3.32 

Aspect  Ratio,  AR 

2.797 

4.277 

2.730 

MAC  (m) 

2.643 

1.180 

1.550 

Incidence,  i(deg) 

1.75 

0 

0 

Reference  Pt:  FS 

626.4 

1120.1 

1079.2 

BL 

195.6 

0 

0 

WL 

280.9 

295.4 

400.6 

Sweep  Angle,  A.  (deg) 

U 

34.0 

32.9 

40.2 

Taper  Ratio,  A 

0.400 

0.201 

0.230 

Body  Diameter  to  Span  Ratio,  d/b 

0.230 

0.160 

0.  163 

2-D  Lift  Curve 

Slope,  aD  (1/deg) 

0.106 

0.109 

0.110 

Zero-lift  Drag 

Coefficient,  Cp^ 

0.0073 

0.0091 

0.0099 

Zero-lift  Moment 

Coefficient,  CL, 

o 

-0.0065 

0 

0 

Max  Lift  Coefficient,  (Ct  ) 

m3X  W(B) 

0.870 

0.980 

1.080 

Angle  of  Attack  at  Max  Cl,  a 

Lmax 

23.0 

24.3 

33.4 

Wing  Downwash:  c/Clw  = 14.07  deg 

Wing  Flap  Effects:  CL0^f  = 0.0058/deg 

ACt  /6f  = 0.0026/deg 
Lmax  1 

ACDo/6f  = 0.00078/deg 

Cm  = -0.0005/deg 
6f 

Wing  Lateral/Directional  Characteristics: 

CL  - 0.060/rad  ACn  /C,2  = 0.055/rad 

B L 
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ac1b/cl  = -0.130/rad 

AClr/cL  * 0.230/rad 

Ci  = -0.211/rad 
P 


TABLE  A-1I1  (Continued) 

AC^/C2  = -0.030/rad 

ACnr/CD  = -0.420/rad 

ACn  /CT  = -0.140/rad 
P L 

ACn  / CD  =18.0 

p ua 


39 , 000  newtons  maximum  installed  thrust.  The  thrust  may  be  vectored  from 
horizontally  aft  (9j  = 0°)  to  five  degrees  forward  of  vertical  (0j  = 95°). 

Nozzle  location  and  assumed  maximum  thrust  per  nozzle  are  listed  in  Table  A-IV. 
The  Pegasus  engine  has  contra-rotating  spools  so  that  its  angular  momentum 
is  negligible.  The  maximum  inlet  mass  flow  was  taken  to  be  181.45  kg/sec  and 
the  resulting  ram  drag  was  assumed  to  act  at  the  centroid  of  the  inlet  area 
(FS  261.2;  WL  246.4). 

Propulsion  induced  aerodynamics  were  approximated  using  the  wind 
tunnel  data  of  reference  (k) . The  basic  data  used  are  shown  in  Figure  A-3 
and  the  functional  approximations  to  these  data  are  given  by  equations  A-l 
through  A-3  where  velocity,  V,  is  assumed  to  be  in  knots 

ALf/T  = -.025  + (.0011944  V - .00005278  V2  + .000000201  V3)  sin  (A-l) 

AD/T  = -.00001386  V + .000008565  V2  (A-2) 

AM/T  = -.15  + (.017082  V + .00003418  V2)  sin  Oj  (A-3) 


Nozz le 
Locat ion 

FS 

(cm) 

BL 

(cm) 

Right  Forward 

465.2 

103.8 

Left  Forward 

465.2 

-103.8 

Right  Aft 

656.9 

79.6 

Left  Aft 

656.9 

-79.6 

TABLE  A-IV 
ENGINE  DATA 


WL 

(cm) 

Toe-out 

(deg) 

Max  Thrust 
(N) 

246.4 

5.0 

14,768 

246.4 

5.0 

14,768 

240.4 

12.5 

14,768 

240.4 

12.5 

14,768 

Aircraft  control  is  accomplished  by  a combination  of  aerodynamic  and 
reaction  controls.  The  aerodynamic  controls  consist  of  conventional  ailerons, 
rudder  and  stabilizer.  The  effectiveness  and  gearing  of  these  surfaces,  as 
modelled,  is  presented  in  Figure  A-4.  Reaction  control  is  achieved  by 
exhausting  engine  bleed  air  through  nozzles  at  the  aircraft  extremities. 

The  locations  of  these  nozzles  are  shown  in  Figure  A-5.  The  pitch  control 
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FIGURE  A-4: 


XV- 6A  AERO  CONTROL  EFF 
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nozzles  (1  and  2)  thrust  only  upward,  the  yaw  control  nozzle  (3)  thrusts  both 
left  and  right  and  the  roll  control  nozzles  (4  and  5)  thrust  both  up  and 
down.  For  the  purposes  of  this  study  the  roll  nozzles  were  assumed  to  thrust 
only  upward  but  with  twice  the  level  of  actual  thrust  so  that  the  moment 
produced  would  be  representative.  Variation  of  RCS  thrust  (at  maximum  engim 
thrust)  with  controller  deflection  is  shown  in  Figure  A-6.  Since  the  RCS  uses 
engine  bleed  air  to  produce  thrust,  nozzle  thrust  will  vary  with  the  main 
engine  thrust.  This  variation  was  approximated  by  equation  A-4. 

%max  RCS  thrust  = 2.393  ( y^  ) - 0.01187  (y—-)2  (A-4) 

where  T is  main  engine  thrust  in  newtons. 

This  completes  the  description  of  the  input  data  used  to  model  the 
XV-6A.  The  model  and  methods  described  in  the  body  of  this  report  were  used 
to  calculate  various  static  and  dynamic  characteristics  of  the  vehicle  for 
comparison  with  available  flight  test  data  taken  from  references  (j),  (k) , 
and  (1).  Results  of  these  comparisons  are  documented  in  the  following. 

Power  off  aerodynamic  characteristics  of  the  model  were  calculated 
both  with  and  without  the  horizontal  stabilizer  over  a range  of  angles  of 
attack  and  sideslips  for  comparison  with  wind  tunnel  results  measured  at 
NASA  Langley  on  a 1/6  scale  model.  Figure  A-7  presents  the  lift  and  moment 
coefficients.  Correlation  is  very  good  for  and  reasonably  good  for  Cm, 

both  tail  off  and  tail  on.  Roll  and  yaw  moment  coefficients  versus  sideslip 
for  two  angles  of  attack  are  presented  in  Figure  A-8.  Again,  correlation  is 
reasonably  good  up  to  sideslip  angles  of  24  degrees.  The  lack  of  measured 
data  precluded  comparisons  at  higher  angles  of  attack  and  sideslip. 

Aircraft  pitch  attitude  and  nozzle  angle  data  were  available  from 
flight  tests  for  longitudinal  trim  from  0 to  240  knots  in  level  flight. 

Figure  A-9  compares  the  estimated  values  with  flight  data.  For  speeds  below 
140  knots  the  aircraft  pitch  attitude  was  constrained  to  be  10  degrees  and 
nozzle  angle  and  longitudinal  control  were  varied  to  obtain  trim.  Above 
140  knots,  the  nozzles  were  fixed  at  0j  = 0 degrees  and  pitch  attitude 

was  allowed  to  vary  for  trim.  The  stabilizer  was  used  for  trim  control 
throughout  the  speed  range  and  was  augmented  by  the  RCS  when  the  nozzle 
angle  was  greater  than  30  degrees.  Correlation  of  nozzle  angle  and  pitch 
attitude  was  very  good.  Due  to  the  unavailability  of  corresponding 
longitudinal  control  position  data  from  flight  test  a comparison  for  this 
parameter  could  not  be  made. 

Estimated  directional  control  effectiveness  is  compared  with  flight 
test  data  in  Figure  A-10.  Engine  nozzle  angle  was  preset  at  the  flight 
test  value  and  the  trim  conditions  were  calculated  for  specified  sideslip 
angles  for  speeds  from  0 to  220  knots.  The  rudder  pedals,  lateral  stick 
and  bank  angle  were  allowed  to  vary  to  achieve  the  specified  trim.  From 
these  data  6/ 6 r was  calculated  and  compared  with  the  flight  test  derived 

values.  Correlation  was  excellent  considering  the  lack  of  data  in  the  model 
for  lateral /directional  propulsion  induced  aerodynamic  effects. 
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FIGURE  A- 8 : XV-6A  LATERAL/DIRECTIONAI. , POWKR-OFF  AERODYNAMICS 
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Hover  control  response  (attitude  and  rate)  was  estimated  for  step 
control  Inputs  about  each  of  the  three  axes  (roll,  pitch,  and  yaw). 

Aircraft  rate  and  attitude  one  second  after  the  initiation  of  the  step 
input  were  calculated  for  comparison  with  similar  quantities  measured  in 
flight.  Figure  A-ll  compares  the  results.  Correlation  is  again  very  good 
for  control  deflections  up  to  2.5  cm.  It  is  interesting  to  note  that 
pitch  response  is  larger  for  forward  (nose  down)  inputs  than  for  aft 
(nose  up)  inputs.  This  is  characteristic  of  the  XV-6A  RCS. 

Pitch  control  response  was  measured  in  flight  and  documented  in 
reference  ( j ) for  forward  flight  speeds  up  to  220  knots.  The  technique 
used  was  to  apply  a slow  sinusoidal  control  doublet  (2  second  period)  to 
the  aircraft  and  measure  the  first  maximum  slope  of  the  pitch  rate  response. 
This  results  in  an  apparent  control  sensitivity  which  includes  the  effect 
of  aerodynamic  pitch  damping.  A similar  approach  was  taken  analytically 
using  the  time  response  of  the  model  to  identical  control  inputs.  The 
results  are  shown  in  Figure  A-12.  Correlation  is  fair  over  the  velocity 
range  considering  the  possible  inaccuracies  in  determining  control  sensitivity 
in  this  manner. 

Aircraft  angular  frequency  response  in  hover  about  all  three  axes  is 
presented  in  Figure  A-13.  Sinusoidal  control  inputs  of  2.5  cm  amplitude 
were  applied  to  the  model  at  frequencies  ranging  from  0.4  to  2.0  Hertz  as 
was  done  in  flight.  The  corresponding  rate  response  for  each  was  calculated 
and  normalized  by  the  amplitude  of  the  control  input.  The  estimated 
frequency  response  correlated  fairly  well  with  that  measured  in  flight. 

A further  examination  of  the  flight  data  indicated  that  the  pitch  response 
data  were  taken  at  a forward  speed  of  10  knots  rather  than  at  zero  knots. 

The  pitch  response  was  estimated  for  the  10  knot  case  and  the  correlation 
improved  as  shown  on  the  figure. 

The  primary  stability  derivatives  were  estimated  for  the  XV-6A  for 
trim  speeds  from  0 to  100  knots.  The  vehicle  was  trimmed  at  8 degrees  angle 
of  attack  in  level  flight  and  the  estimated  derivatives  were  normalized 
by  mass  or  the  moment  of  inertia  about  the  corresponding  axis.  The 
estimated  results  are  compared  with  National  Research  Council  of  Canada 
estimates  for  the  P1127  (reference  (1)),  McDonnell  Aircraft  Corporation 
estimates  for  the  AV-8A  and  limited  results  derived  from  flight  tests  of 
XV-6A.  Figures  A-14  and  A-15  present  the  resulting  data  comparisons. 

Longitudinally,  the  estimates  appear  reasonable  for  all  derivatives 
except  and  Zu.  Laterally,  the  estimates  compare  favorably  with  the 

AV-8A  data  for  all  derivatives  except  Lv,  Lr,  and  Np.  The  NRC  estimated 

derivatives  do  not  compare  with  either  the  AV-8A  data  or  the  model  estimates. 
It  is  difficult  to  draw  firm  conclusions  on  the  accuracy  of  the  subject 
method  due  to  scatter  in  the  available  data  which,  for  the  most  part,  are 
also  estimates  and  are  not  substantiated  by  flight  test  data. 

In  general  , the  method  described  in  the  body  of  this  report  appears  to 
yield  reasonably  accurate  results  compared  to  available  flight  test  data — 
at  least  for  Jet-lift  vehicles  of  the  Kestrel  type.  Further  validation  is 
necessary  to  extend  the  credibility  level  of  the  method  to  other  classes 
of  V/STOL  vehicles. 
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FIGURE  A- 12 : XV-6A  PITCH  CONTROL 
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